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Recent radio observations support a picture for star formation where there is accretion of matter onto a central 
protostar with the ejection of molecular outflows that can affect the surrounding medium. The impact of a supersonic 
outflow on the ambient gas can produce a strong shock that could accelerate particles up to relativistic energies. A 
strong evidence of this has been the detection of non-thermal radio emission coming from the jet termination region 
of some young massive stars. In the present contribution, we study the possible high-energy emission due to the 
interaction of relativistic particles, electrons and protons, with the magnetic, photon and matter fields inside a giant 
molecular cloud. Electrons lose energy via relativistic Bremsstrahlung, synchrotron radiation and inverse Compton 
interactions, and protons cool mainly through inelastic collisions with atoms in the cloud. We conclude that some 
massive young stellar objects might be detectable at gamma-rays by next generation instruments, both satellite-borne 
and ground based. 
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1. Introduction 

The mechanism of formation of massive stars (M > 
8Mq) remains one of the open questions in astro- 
physics. It is known that these stars originate inside 
giant and massive molecular clouds but the sequence 
of processes that take place during the formation of 
the star are mostly unknown. It has been suggested, 
for example, that the coalescence of various proto- 
stars in the same cloud can lead to the emergence of 
a massive star [1]. Massive stars appear in massive 
stellar associations where cloud fragmentation seems 
to be common. Alternatively, a massive star could 
form by the collapse of the core of a molecular cloud, 
with associated episodes of mass accretion and ejec- 
tion, as observed in low-mass stars [2] . In such a case, 
the effects of jets propagating through the medium 
that surrounds the protostar should be detectable. 

Until now, the formation of stars has been mostly 
associated with thermal radio and X-ray emission. 
However, non-thermal radio emission has been de- 
tected in some massive star forming regions. This is 



a clear evidence that efficient particle acceleration is 
occurring there, which may have as well a radiative 
outcome at energies much higher than radio ones. 

In the present contribution, and based on recent 
multiwavelength observations, as well as reasonable 
physical assumptions, we show that massive proto- 
stars could produce a significant amount of radiation 
in the gamma-ray domain, because of the dense and 
rich medium in which they are formed. 

2. Non-thermal radio sources 

In the last years, synchrotron radiation have been ob- 
served from some regions where massive stars form. 
This emission is associated with outflows emanat- 
ing from a central protostar. In what follows, we 
briefly describe some of these non-thermal sources 
that could be potential emitters of gamma-rays. 

2.1. IRAS 16547-4247 

The triple radio source associated with the protostar 
IRAS 16547-4247 is one of the best candidates to 
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produce gamma-rays. This system is located within 
a very dense region (i.e. densities n « 5 x 10 5 cm~ 3 ) 
of a giant molecular cloud located at a distance 
of 2.9 kpc. The luminosity of the IRAS source is 
L = 6.2 x lO 4 ^©, possibly being the most luminous 
known YSO associated with collimated thermal jets. 

ATCA and VLA observations [3, 4] have shown 
that the southern lobe of this system, of size ~ 
10 16 cm, has a clear non-thermal spectrum, with an 
index a ~ 0.6 (S u oc v~ a ). The specific flux of this 
source is 2 mJy at 14.9 GHz and the estimated mag- 
netic field is B ~ 2 x 1(T 3 G [5]. 

2.2. Serpens 

The Serpens molecular cloud is located at a distance 
of ~ 300 pc. One of the two central dense cores of 
this cloud is a triple radio source, composed by a 
central protostar (IRAS 18273—0113) and two lobes. 
The northwest (NW) hot-spot is connected with the 
central source by a highly collimated thermal jet, 
whereas the southeasth (SE) is separated and bro- 
ken into two clumps. The luminosity of the source 
IRAS 18273-0113 is L ~ 3OOi and the parti- 
cle density at the center of the molecular cloud is 
n ~ 10 5 cm~ 3 . 

The observed radio emission [6, 7] detected from 
the central and NW sources has a spectral index 
a w —0.15 and a « 0.05, respectively. This emission, 
of a luminosity ^2 — 3 mJy, is produced via thermal 
Bremsstrahlung. However, the radiation produced in 
the SE lobe seems to be non-thermal (a = 0.3), likely 
produced via synchrotron emission. The specific flux 
of this source is 2—5 mJy. The equipartition magnetic 
field estimated in the SE lobe is -B eq uip ~ 10~ 3 G [6]. 

2.3. HH 80-81 

The famous Herbig-Haro objects called HH 80—81 
are the south component of a system of radio sources, 
located in the Sagitarius region, at a distance of 
1.7 kpc. The central source has been identified 
with the luminous (L = 1.7 x 10 4 L Q ) protostar 
IRAS 18162-2048. HH 80 North is the northern 
counterpart of HH 80—81. The velocity of the jet has 
been estimated as v ~ 700 km s _1 , allowing to derive 
a dynamical age for the system similar to 4000 yr. 

Radio observations carried out with the VLA in- 
strument [8] showed that the central source has a 
spectral index a ~ 0.1, typical of free-free emission, 



whereas HH 80-81 and HH 80 North are likely non- 
thermal sources, with spectral index a ~ 0.3. The 
specific flux measured at a frecuency of 5 GHz is 
F v ~ 1 — 2 mJy and F v w 2.4 mJy for the sources 
HH 80-81 and HH 80 North, respectively. At this 
frecuency, the angular size of the north and the south 
components are ~ 6". 

In addition, HH 80—81 system is a powerful 
source of thermal X-rays with a luminosity Lx ~ 
4.3 x 10 31 erg s _1 [9]. 

2.4. W3(OH) 

Other interesting source to study is the system com- 
posed by an H2O maser complex and the Turner- 
Welch (TW) source in the W3 region [10]. The cen- 
tral source of this system is a luminous (L ~ 1O 5 L ) 
YSO and the mean density of cool particles is n ~ 
4 x 10 4 cm- 3 . The distance to W3(OH) is 2.2 kpc. 

Continuum radio observations [10] show the 
presence of a sinuous double-sided jet, emanating 
from the TW source. The observed radio flux, from 
1.6 to 15 GHz, is in the range 2.5-0.75 mJy, and the 
spectral index of the observed emission is clearly non- 
thermal: a = 0.6. The inhomogeneous synchrotron 
model proposed by Reid et al. in Rcf. [11] pre- 
dicts for the emitting jet a density of relativistic 
electrons n e (j,r) = 0.0687~ 2 (r/r )~ 1 ' 6 cm~ 3 and 
a magnetic field B{^,r) = 0.01(r/r )~ ' 8 G, where 
ro = 6.6 x 10 15 cm, and r and 7 are the distance to 
the jet origin and the electron Lorentz factor, respec- 
tively. Unlike in previous cases, here the non-thermal 
radio emission comes from the jet and not from its 
termination region. Non-thermal jets associated to a 
YSO are uncommon. 

3. Acceleration of particles and losses 

The non-thermal radio emission observed in some 
massive star forming regions is interpreted as syn- 
chrotron radiation produced by the interaction of 
relativistic electrons with the magnetic field present 
in the cloud, being typically i? c ioud ~ 10~ 3 G [12]. 
These particles could be accelerated at a shock, 
formed in the point where the jet terminates, via 
diffusive shock acceleration [13]. The acceleration ef- 
ficiency, characterized by 77, is related to the veloc- 
ity of the shock. Using the values of the velocities 
given for the sources described in the previous sec- 
tion, and assuming Bohm diffusion, values for 77 of 
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~ 10~ 6 - 10~ 5 are obtained. 

Particles accelerated up to relativistic energies 
interact with the different fields present in the 
medium. As noted at the beginning of this section, 
electrons radiate synchrotron emission under the am- 
bient magnetic field B. In addition, particles, elec- 
trons and protons, can also interact with the cold 
matter in the jet termination region (via relativistic 
Bremsstrahlung the leptons, and proton- nuclei col- 
lisions the protons). In addition, electrons interact 
with the background field of IR photons of the pro- 
tostar, of energy density w pri , via inverse Compton 
(IC) scattering. 

Using the following parameter values: n = 5 x 
10 5 cm- 3 ; B = 2.5 x 1(T 3 G; and u ph = 3.2 x 
1CT 9 erg cm" 3 given for IRAS 16547-4247, [3, 5], 
we estimate the cooling time of the main leptonic 
processes in this scenario. As seen in Fig. 1, rela- 
tivistic Bremsstrahlung losses are dominant up to 
~ 10 GeV. In addition, it is possible to see from this 
figure that the maximum energy achieved by elec- 
trons is £;™ ax ~ 1 TeV and is determined by syn- 
chrotron losses, being IC losses negligible. 
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Fig. 1. Energy loss and acceleration rates for electrons in the 
IRAS 16547-4247 southern lobe. 

Protons can be accelerated by the shock in the 
same way as electrons and interact with cold particles 
present in the cloud. The maximum energy achieved 
by protons is higher: E™ ax ~ 10 TeV. In pp inter- 
actions, besides 7-rays, secondary electron-positron 
pairs are produced. These secondary particles will ra- 
diate by the same mechanisms as primary electrons 
(i.e. synchrotron radiation, IC scattering and rela- 
tivistic Bremsstrahlung) . 



4. High-energy emission 

In order to calculate the non-thermal spectral en- 
ergy distribution (SED) of a massive protostar, the 
magnetic field of the cloud and the distributions of 
relativistic particles, n(E), arc needed. To obtain the 
magnetic field B and the normalization constants of 
the particle distributions, we use the standard equa- 
tions given by Ginzburg & Syrovatskii [14] for the 
observed synchrotron flux and assuming equipatition 
between the magnetic and the relativistic particle 
(primary and secondary pairs and protons) energy 
densities: 

B 2 

— = u e + u p + u e ±, (1) 

In Eq. (1), the following relationships are implicit: 
Up = au c and u e ± = f u p . In the first relation a is 
a free parameter (here fixed to 100, as it is the case 
for cosmic rays) and in the second relation / can be 
estimated using the average ratio of the number of 
secondary pairs to 7r°-decay photons in Ref. [15]. In 
Fig. 2, we show the computed broadband SED for 
the parameters of the source IRAS 16547-4247. 

5. Discussion 

In this work we show that, if the source is located at 
few kpc, the high-energy emission may be detected 
by GLAST and even by forthcoming Cherenkov tele- 
scope arrays after long enough exposure. This opens 
a new window to the study of star formation and 
related processes. Also, determinations of the parti- 
cle spectrum and its high-energy for different sources 
with a variety of environmental conditions can shed 
light on the properties of galactic, supersonic out- 
flows, and on the particle acceleration processes oc- 
curring at their termination points. 

Radio observations already demonstrate that rel- 
ativistic electrons are produced in some sources. Ac- 
cording to the presence of non-thermal emission de- 
tected at cm- wavelengths and IR observations of the 
protostars emission we can suggest several good can- 
didates to be targeted by GLAST. These objects are 
IRAS 16547-4247 [5], HH 80-81 [8], W3 [10, 11] and 
the multiple radio source in Serpens [6, 7] 

To conclude, we emphasize that massive YSO 
with bipolar outflows and non-thermal radio emis- 
sion can form a new population of gamma-ray 
sources that could be unveiled by the next gener- 
ation of 7-ray instrumentation. 
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Fig. 2. Spectral energy distribution for the south lobe of the YSO embedded in the source IRAS 16547—4247. The radiative 
contribution of secondary pairs (2) produced via Tr^-decay is shown along with the contribution of primary electrons (1). 
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